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Abstract: We investigate how nontrivial topology affects the entanglement dynamics 
between a detector and a quantum held and between two detectors mediated by a quantum 
held. Nontrivial topology refers to both that of the base space and that of the bundle. Using 
a derivative-coupling Unruh-DeWitt-like detector model interacting with a quantum scalar 
held in an Einstein cylinder (space) xRi (time), we see the beating behaviors in the 
dynamics of the detector-held entanglement and the detector-detector entanglement, which 
distinguish from the results in the non-compact (H-l) dimensional Minkowski space. The 
beat patterns of entanglement dynamics in a normal and a twisted held with the same 
parameter values are different because of the difference in the spectrum of the held modes. 
In terms of the kinetic momentum of the detectors, we hnd that the contribution by the 
zero mode in a normal held to entanglement dynamics has no qualitative difference from 
those by the nonzero modes. 
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1 Introduction 

The last decade has seen a rapid increase of research in relativistic quantum information 
(RQI) (see, e.g., [1]), which studies the relativistic features of quantum information with 
“relativistic” referring to the localized quantum objects in relativistic motion, the relativistic 
nature of quantum fields and the properties of spacetime structures. In particular, the 
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third aspect concerning spacetime properties ranges from the rudimentary yet necessary, 
such as time-dilation and hypersurface-slicing considerations [2, 3] at the level of special 
relativity, to spacetime curvature and topology effects [4, 5] often considered in the context 
of general relativity. The present work will focus on how nontrivial topology affects the 
entanglement dynamics of a system of (= 1 or 2 here) detectors, which are localized atom¬ 
like quantum objects with internal degrees of freedom sensitive to the field they interact 
with. Here, nontrivial topology refers to the base space, turning a flat, noncompact two 
dimensional Minkowski spacetime R} = R^ (space) xRi (time) into a spatially-compact 
(space) xRi (time), and also the bundle, whose fibre is the range of the field and a field 
configuration is a cross section - we are interested in seeing how entanglement dynamics 
for a twisted field is different from a normal (untwisted) field. 

Before proceeding further we hasten to point out that although the mathematical na¬ 
ture of these investigations is obvious, they are physically relevant to actual experimental 
setups: boundary effects pertain to the effects of a mirror or a dielectric slab [6] on en¬ 
tanglement, while a spatial topology refers to a toroidal cavity [7], both are obviously 
essential components in quantum experiments. On the theoretical side, quantum field the¬ 
ory in spacetimes with boundaries or with nontrivial topology was investigated in the late 
70s by Isham, Dowker, DeWitt, and co-workers [8-10]. The prototype spacetime x Ri 
was referred to as the Einstein cylinder (presumably inspired by the “Einstein Universe” 
with topology x Ri, a static closed Robertson-Walker Universe). The new aspect in our 
present investigation is entanglement dynamics in this topologically nontrivial spacetime 
[11], for both untwisted and twisted quantum fields 

Let us ponder upon how the overall features of entanglement dynamics between detec¬ 
tors and fields or between two detectors in a common field would be modified in a circle 
(e.g., a narrow micro-toroidal cavity) as compared to a straight line, as investigated before 
in e.g., [14]. 

Topology of the base space In R^ space, dynamics of the detector-field or detector- 
detector entanglement shows no periodicity at large time scales, we expect entanglement 
would just grow or decay and saturate. In S^, the space is finite and compact, so the 
retarded field emitted by the detector will after one period around the circle return to affect 
the detector and do so periodically. In a fully enclosed space there is no true dissipation, 
only apparent one on a time scale short compared to the recurrence time, as well illustrated 
in the Kac ring [15], a ID closed harmonic chain, but the energy apparently “lost” to the 
rest of the chain will be replenished after one period This of course is the advantage for 
information processing using micro-toroidal cavities (the space in realistic situations is of 
course not fully enclosed, as there will be input and output laser activities). For quantum 
systems the time lag between successive cycles generates interference which shows up as 
beats. Unlike in R^, where equilibrium or steady state can be found to exist at late times, 
in S^, beating in entanglement distinguishes its behavior. 

^The effect of the topology of the state-vector space on the decoherence and entanglement processes has 
been investigated for some time. See, e.g., [12, 13] 

^Entanglement propagation in a quantum Kac ring has been studied by [16]. 
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Topology of the bundle Here the nontrivial topology exists in twisted fields, thus the 
comparison we aim at is between the entanglement dynamics of normal (untwisted) and 
twisted fields. In S^, the untwisted and twisted fields possesses different eigen-modes, and 
hence different beat patterns. The normal field possesses a zero mode whereas for twisted 
fields there is no zero-mode. Note that the existence of the zero mode is a topological, not 
a geometric effect, in the sense that for normal fields different sizes of as in a toroidal 
cavity will all possess a zero mode. Likewise there is no zero mode for twisted fields for all 
sizes of cavity. 

We mention two related work [4, 11] which uses the time-dependent perturbation theory 
with a switching function type of interaction to study the finite-time response of a single 
2-level Unruh-DeWitt detector and the quantum entanglement of two such detectors [11] 
in locally flat (3-1-1) dimensional spacetimes with different topologies. The lowest-order 
results do depend on the field spectrum, which in turn depend on the spatial topology, 
though the differences in different topologies are small. However, their lowest-order results 
do not show beating behavior. The main reason is that the influences by the echoes of 
a detector itself, as well as the retarded influences by the other detectors separated at a 
distance, arise from the higher order terms in the coupling expansion. 

Below, we use a nonperturbative method to study the nonequilibrium dynamics of 
a harmonic-oscillator detector-field system with all the mutual influences included in our 
consideration, in order to capture the full extent of the topological effect in entanglement 
dynamics. The setup in our study is introduced in Section 2. We consider a Unruh- 
DeWitt-like (UD’) detector theory with derivative-coupling [17-22] because a usual detector 
minimally coupled to a scalar field in (H-l) dimensional Minkowski space has undesirable 
non-positive-definite energy and higher-derivative radiation-reaction. ^ In Section 3 we 
consider one detector interacting with a quantum field in an Einstein cylinder and compare 
its entanglement behavior with earlier results obtained for Minkowski space R^. In Section 
4 we consider the entanglement dynamics of two identical UD’ detectors in the presence 
of a quantum field in the same x Ri spacetime for both normal and twisted fields. 
Our findings are summarized in Section 5. More details about the two-point correlators 
of a UD’ detector in R^ are presented in Appendix A. In our numerical results, we set 
c = h = 1. 

2 Model 

Consider a (1+1) dimensional flat spacetime with topology in space and Ri in time, 
namely, the Einstein cylinder. The metric is given by 

ds'^ = —dt^ + dx‘^, (2.1) 

where x = R<f with radius R a positive real constant and the azimuthal angle ip G (—tt, tt]. 
The circumference of the circle is thus L = 2ttR. We will also refer to the extended space 

® An additional advantage when using time-dependent perturbation theory is that the infrared behavior 
of the response function of a detector with derivative-coupling is better than that with minimal-coupling 
in (U-l) dimensions, see [23]. 
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x,x' € obtained by identifying the points x to x mod L (see Figure 1). Consider placing 
a finite number of the derivative-coupling Unruh-DeWitt-like (UD’) [17-20, 22] detectors 
with the internal harmonic oscillators Qd coupled to a common massless scalar field in 
the above spacetime, described by the action 




where the interaction action is 


^ J d‘^xda<^d^^ + Z] ^ [(^dQd)^ - ^iQd] 


+ St 


5’/ = -^ A J dr^Qddd J d‘^x^{t,x)5'^{x°‘ - 


( 2 . 2 ) 


(2.3) 


with x“ = (t,x), a = 0,1, the detector label d = ^ for one-detector case, d = A, B for 
two-detector case, and = d/dr^. 2 ;§(Tjj) are the prescribed trajectory of the detector d. 
The canonical momenta conjugate to and are 

6S 


PM = 




= ddQd i^d) > 


(2.4) 


n(t,x) = = dt^{t,x) - Z^/ dTdQdi'^d)Md)^^ix'" - 4(^d))> (2-5) 

respectively, where Then, after a Legendre transformation, one can write down 

the Hamiltonian as 

= Z 2v^{t) + ‘^dQd('’'d(^))} + 

-\ 2 




n^W + Z^ / drdQd{Td)vd{^d)^^{x'" - Zd{Td)) 

J 'J 


+ > , ( 2 . 6 ) 


which is parametrized by the time variable x*^ = t of the observer’s frame and defined on 
the whole time-slice x G associated with t. 

Alternatively, one may adopt the interaction action 


S'i = '^X J dXdddQd J d'^x<^{t,x)5‘^{x°^ - Zd{Td)), 


(2.7) 


plus surface terms evaluated at the initial and hnal moments from (2.3), which gives the 
same Euler-Lagrange equations for the dynamical variables and ^x- Starting with (2.7), 
the canonical momenta conjugate to and 4? become 

Pdi^d) = ^dQd (''“d) + A‘h (2;d(Td)) , (2.8) 

U'{t,x) = dt^{t,x), (2.9) 

respectively, and so the Hamiltonian reads 


x(t} = E 


1 2 


Pdi^) ~ (^) + ‘^dQd('''d(^)) f + 


d 2^11(t) 

1 [ <ix{nJ(«) + [ 8 A(*)|"}. 


( 2 . 10 ) 
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Figure 1. (Left) The setup. The harmonic oscillators Qa and Qb are the internal degrees of 
freedom of detectors A and B located at a; = 0 and a: = L/2 = Ett, respectively. At f = 0 we let 
the UD’ detectors begin to couple with the field. The orange and blue solid lines are the worldlines 
of the two detectors, and the dashed lines represent the fronts of the retarded influences sourced 
by the detectors after the interaction is switched on. (Middle) An equivalent setup in the extended 
coordinates in space by identifying the points a; to a; mod L . (Right) Similar diagram for the 
cases with a single detector Qa only. 


While the values of the Hamiltonian H' and the momentum are invariant under a 

shift —)• <ha;+C with a constant C, such a global symmetry is enough to make the value of 
Pj = Qd + ('^d('''d)) “ga-uge” invariant. Since we are looking at the reduced state of 

the detectors in this paper, we prefer to work with (2.3) to get rid of some weird dynamical 
behaviors of P'^ which could be “gauged away” by the surface terms mentioned below (2.7). 
Later we will see that the two point correlators {P'^) and the corresponding uncertainty 
functions will have indefinite growths in time, while (Pj) of the kinetic momentum is 
well-behaved. 

Suppose the combined system of the detectors and the field is initially in a Gaussian 
state, which could be pure or mixed. Then the quantum state will always be Gaussian in 
its history because of the linearity of the interaction. A Gaussian state is fully determined 
by the two-point correlators of the dynamical variables of the combined system, which can 
be obtained by taking the expectation values of two operators, evolved in the Heisenberg 
picture, with respect to the initial state (see, e.g., [14]). 

The Heisenberg equations of motion for the operators read 

^dQd('^d) +‘^dQd('''d) = “-^^d^’(2^d('''d))) (2-11) 

-□6(x) = [ drdddQdi-rd)^^ - Zdi^d)) ^ (2-12) 

where □ = —d^ + 5^. Suppose the detector-field coupling is suddenly switched on at the 
initial moment t = 0 (when = 0 for all detectors). By virtue of the linearity of the 
system, the operators will evolve into linear combinations of the operators initially defined 
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at t = 0 as 



9d (^d)ad' +9d*('^d)a 


d'* / 


7^ Q'd('''d)^fc + 9d*('''d)^I ) (2-13) 


= 



2CcJd - 



+ ^x*i't)^d' + y~! \l 77?r ^x{'t)ik + 4’x*{tWk )(2-14) 


with tdfc = WfcL = |/c|L for /c / 0 and Wfc = 1 for A: = 0. Here q^', and (j)^ are the 

(c-number) detector-detector, detector-field, field-detector, and field-field mode functions, 
respectively, and a|j are the lowering and raising operators for the free detector d, while 
bf^ and are the annihilation and creation operators for the free field mode with wave 
number k. For the normal (untwisted) field, we take k = kn ^ n/R = 27rn/L, n G Z 
to satisfy the periodic boundary condition and for the twisted field, 

we take k = k'^ = [n — (l/2)]/i?, n G Z to satisfy the anti-periodic boundary condition 
(analogous to the Mobius band) [8, 10]. There is one zero mode 
(/cq = 0 and Ukf) = 0 when n = 0) for the untwisted held, and no zero mode for the twisted 
held. For the zero mode, we dehne bo = {^ko + with the initial zero-mode held 

operator and its conjugate momentum operator satisfying the equal-time commutation 
relation [$fcQ,nfco] = ih. 

Since the detectors and the held are free before the detector-held coupling is switched 
on at t = Td = 0, the initial conditions for the mode functions are set to be <?if](0) = 
Mxi^) = for /c / 0, ^(^'(O) = 5^', 545;^'(0) = -iw^^d > and 0^(0) = (9t^^(0) = 

9d(^) ~ ^tQdi^) ~ zero-mode component of the free held [24], one has ^A:o(0 = 

, so the initial conditions 

are set to be <?i>^°(0) = 1 and d-i.(j)^°{0) = —ijL. The momentum operators conjugate to 
and at t > 0 can be obtained straightforwardly by inserting (2.13) and (2.14) to the 
operator version of (2.4) and (2.5). 

In terms of the above expansion, Eq. (2.12) yields 


l>fc„(0) + {t/L)Uko{0) = vW2 (1 - + (1 + {i/mb 


- = [ dmdQdim^ - 2d(^d)) (2-15) 

d d 

where /i = d,k, and the mode functions with different /i’s are decoupled. Eq. (2.15) 
implies c/)x{t) = 4> x^{t) + 4> x^{t), where the homogeneous solutions are = 0, = 

^-iuj^.t+tkx A: ^ 0, and (/>^°'^°(t) = 1 — {i/L)t for the zero mode, and the inhomogeneous 
solutions read 

fll poo 

= ^^dG'ret(a:“;4(rd))(9dgd(^d) (2-16) 

d do 

with the retarded Green’s function Gret of the held. When all of the held modes are 
considered, the retarded Green’s function Gret for the massless scalar held can be written 
as 

Gret(A, x; t' , x') = - e^6 \t + x — {t' + x' + nL)] 6 \t — x — {t' — x' — nL)] , (2.17) 

^nez 
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for X, x' G (—L/2, L/2]. We take e = 1 for the untwisted field, and e = —1 for the twisted 
field. In either case, the retarded Green’s function looks the same as the one in Minkowski 
space around the source point (f,x) = {t',x'). One can verify that —OGret{t,x‘,t',x') = 
Yln&z — t')5{x — x' — nL) in the extended coordinates with x, x' G and identifying 
the points x to x mod L (see Figure 1). Now can be interpreted as vacuum 

fluctuations of the field state, while (j)x ^{t) behave like the retarded helds sourced by the 
point-like detectors. Note that, to be consistent with the expressions (2.17), later we have 
to include the contributions by enough number of the the field modes, namely, the UV 
cutoff cannot be too small. 

For the untwisted field, Eq. (2.17) includes the contribution by the zero mode. If we 
exclude the zero mode entirely, the retarded Green’s function for the untwisted field will 
become [24] 

G'^l{t-,x\t',x) 

= hm^ - t')Im {in + In | 

= - 9 + X — {t' + x' + nL)] 9 — X — {t' — x' — nL)] — ^ 9{t — t'). (2.18) 

nez 

Compared with —□Gj.g^, here —has an extra term —□ [—(t — t')9{t — t')/L\ = —(2/L) 
8{t — t'). This is because the field modes cannot form a complete basis without the zero 
mode. Moreover, is nonzero when \t — t'\ < jx — x'j, which implies that two spacelike 
separated events located at (t, x) and {t', x') can have superluminal signaling, which violates 
causality, even at the classical level. 

3 One-detector case 

Let us start with the simplest case with only one single detector A with natural frequency 
= rio) located at x = 0. To study the influence of the quantum field in the x Ri 
spacetime on the detector we look at the reduced dynamics of A by tracing out the field. 
For this purpose we need to calculate the two-point correlators of the detector. 

Inserting the solutions of (px with (2.16) into (2.11), one obtains the equations of motion 
for the mode functions 

OO 

{d^ + + L?^) q'^it) = “ n'L)d^q^^{t - n'L) (3.1) 

n'=l 

^ f for t < L .g 2 ) 

\ e (9| - 2-idt + Sd^) q^j^{t - L) for t > L, 

with fi = A, kn (untwisted field) or (twisted field), n G Z, and the coupling strength 
7 = A^/4. When t > L, the above expressions for different /r’s have the same appearance, 
while the solutions at 0 < t < L before the first echo hits the detector are generated by 
different driving forces oc dt4>Q^{t). Note that the above equations for the twisted field on 
are equivalent to those for a UD’ detector in the same scalar field and located at the 
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center of an ID cavity made of two perfectly-conducting mirrors separated at a distance of 
L, while the ones for the untwisted field are equivalent to those in a cavity with infinitely- 
permeable mirrors [6], or in a toroidal cavity [7]. 

Although the left-hand side (LHS) of (3.1) describes the continuous evolutions of the 
mode functions, on the right-hand side (RHS) of the equation the influences by the echoes 
come in a discrete fashion. The analytical solutions for these delayed differential equa¬ 
tions can be obtained order by order from those at very early times in principle, but the 
lengths of the expressions for the solutions grow rapidly as the order of the included echoes 
increase. Very soon the analytical solutions will get too complicated to read off any use¬ 
ful information. In this case numerical computation would come in handy and give more 
transparent results. 

Below, the delayed differential equations (3.2) will be solved with the proper initial 
conditions. For t < L, the solutions for (3.2) have the closed form 


^Ait)\o<t<L 




[(D - (Do + fy)) + (D + (Dq + iy)) e”*™] 


0<i<L 


—A / drD ^K{t — T)d^(f) q^{t) 

Jo 

Xlo ( 

2D i y — i(a;-I- D) y — i(u; — D) 


_ g(-7-fir2)i 


A“(*) 


— —A / drQ ^K(t — T)dfd>^Q^(r) 


0<t<L 


iX 

L 


1 — e ^cos sin Dt^ 


(3.3) 


(3.4) 


(3.5) 


with oj = \k\ {k can be kn or k!^), D = y^Dq — y^, and K{x) = e '>'^sinDa:. The above 
early-time solutions for (?^(t) and are identical to those in Minkowski space [25]. 


3.1 Eigen-frequencies 

3.1.1 Untwisted field 

Assume at late times {t ^ L, 1/y) the detector would evolve into a stationary state, when 
the mode function could be written as 

qi{t) « j duqi{u)e"'^\ (3.6) 

For e = 1, inserting the above ansatz into (3.1) or (3.2) with ^ = A yields 

— -|-Dq = —2ya; cot(a;L/2). (3.7) 

The solutions of w are the eigen-frequencies. One can immediately see that |a;| = Dq are 
solutions but only when i? = (1 -|- 2n)/(2Do) (i.e. L = 7r(l -|- 2n)/Do), n = 0,1, 2,3,.... 
More general solutions to (3.7) can be obtained by numerical methods. In Figure 2 (upper- 
right) one can see that when Dq 27r/L, for those |a;| < Dq (|w| > Dq), one has |a;| < n/R 
(Iwl > n/R), n = 1, 2, 3,. 
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When rio ~ n/R with some positive integer n, while in (3.4) is on resonance at 
early times (which, together with other modes, produces some resonant oscillation on top of 
the decaying behavior of the two-point correlators and functions of them, see Appendix A), 
at a larger time scale the detector mode with Qq and the field mode with LOk„ = \kn\ = n/R 
will mix together and generate two dominant eigen-modes ^ with frequencies w ~ i A, 
where 

A RS \/^, (3.8) 

in the weak coupling limit. These two eigen-modes will dominate the late-time behavior 
of q^. In particular, the frequency difference produces the beat of at frequency 2A, 
and so 27r/(2A) ~ is the largest significant time scale in the evolution of the single 

detector system with Qq ^ n/R in the weak coupling limit. This is very different from the 
time scale 7 “^ in the detector theories in Minkowski space. 

Even if the radius R or the circumference L of is larger than the time scale I/ 7 , 
which is the largest time scale of the detector in (H-l) dimensional Minkowski space R}, 
only the early-time {0 < t < L) behavior of the mode function will be similar to those in 
RJ. Once the echoes start to affect the detector {t > L), the higher-order corrections from 
the echoes will be the same order of magnitude as the zeroth order solution. At large time 
scales one can also see the beat feature when Hq ~ n/R for some integer n. 

When IIo is not close to any n/R, n ^ Z, the beating behavior may be suppressed. In 
the example shown in Figure 3, when Oq ~ [n — {l/2)]/R for some n, the only dominant 
eigen-mode has a frequency very close to Hq, and the largest time scale in the evolution 
of the mode function is about L, which is the period of the massless field traveling over 
the space S^. However, if there are more than two dominant modes, the small differences 
between the frequency-differences of the dominant modes may produce beats at an even 
larger time scale (e.g.. Figure 7). 

For the detector-field mode function q^ , inserting a similar ansatz 

q^/{t)^ J didq'‘/{uj)e^‘^\ (3.9) 

into (3.1) or (3.2) with ^ = kn yields 

— -|- Hq = — 27 CU cot(a;L/2) — — w). (3.10) 

The solutions of the eigen-frequencies for q^J/ are exactly the same as those for q^ with the 
same values of the parameters Hq, A, and L. This is evident when comparing Figure 4 with 
Figure 2, and Figure 5 with Figure 3. One can clearly see the beats in Figure 5 while there 
is no significant beat in Figure 3. This is because there are two dominant eigen-modes 
(with frequencies a; ~ Oq and |A:„|) for q^ in Figure 5, rather than one for q^ in Figure 3. 

^In this paper, our “eigen-modes” refer to those in the stationary ansatz for a mode function, e.g., 
(3.6) or (3.9), rather than the mode fnnction itself. We call the eigen-modes “dominant” if the amplitude 
|g(cu)| in the ansatz has the maximum values at the corresponding eigen-frequencies. 
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{y=0.1,n=2.3,R=1.3) 

1^1 jY=0.1, £2=2.3, R= 1.31 



Figure 2. Emergence of the frequency spectrum, after a Fourier transform of from t = 0 to 
t = T ^ L (lower-right). In the upper-right plot, the blue and purple curves represent the LHS and 
RHS of Eq.(3.7) as functions of ujR, respectively. Here the LHS = 0 at the free detector’s natural 
frequency Hq = ±2.9928/i? « 3/R, where a free-held mode is located. These two frequencies mix 
together and generate two new eigen-frequencies ojR « 3.1901 and 2.8089, which are the closest 
two among all of the eigen-frequencies in the frequency space. In the lower-right plot one can 
see that these two eigen-frequencies dominate the long-time behavior of and give the beat of 
in the left plots. The period of the beat can be estimated from the difference of these two 
eigen-frequencies: Tbeat ~ 27r/(2A) = 27r/[(3.1901 — 2.8089)/i?] « 21.4297 (the gray dashed curve 
in the upper-left plot represents (1 -I- cos 27rt/Tbeat)/2 « (1 -I- cos 2At) /2). Compare the upper- and 
lower-left plots one can see that, while the frequency spectrum gets sharper and sharper as the 
duration of interaction T increases, the mode function have similar behaviors from early times all 
the way to late times, if observed in the same small time scale. 

3.1.2 Twisted field 

For e = —1, inserting ansatz similar to (3.6) or (3.9) into (3.1) or (3.2) yields 

—-b Oq = 27 a; tan(t(jL/ 2 ), (3-11) 

—-b Oq = 27a;tan(ti;L/2) — X2Trd{uik'^ — uj), (3-12) 

for ^ = A and /r = k'^, respectively. Again the eigen-frequencies of q^{t) and q^{t) obtained 
from the above two equations are the same. Since tan(a;L/2) = — cot{ujL/2 — irI2f)^ for a; 
large enough and not very close to Qq, the solutions of (3.11) and (3.12) are roughly those 
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Figure 3. When 7 is small and fig is about [n — (l/2)]/i? (n G Z), which is not close to any n/Ji, 
no significant beat can be observed. Here 7 = 0.01, fl = 2.3, and R = 2. The most significant 
eigen-mode has the frequency closest to Hg. 


|T=0.1.n=2.3,R=1.31 (7=0.1,n=2.3, R=1.31 ( v=0.1. £i=2.3. R=1.31 ( Y=0.1. n=2.3, R=1.3 ) 

\qi\^ l94p I9 a|’ ki’l’ 



Figure 4. q^{t) (upper row) and their frequency spectra (lower row), n = 2,3,4, 5 (from 

left to right), with the same values of the parameters as those in Figure 2 (7 = 0.1, H = 2.3, 
i? = 1.3). One can see that the eigen-frequencies are exactly the same as those for the there. 


solutions for (3.7) shifted by tt/L = 1/(27?). (Recall that the modes for the free twisted 
field have = [n — (l/2)]/72.) 

The difference between the spectra of the free fields can alter the beating behavior 
drastically in different fields while the values of the parameters Oq, A, and L are the same. 
For example, for Oq ~ n/R, n G Z, the beating behavior of is not significant in the 
twisted field, in contrast to the clear beats in the cases with the same parameter values 
but in the untwisted field, as those in Figure 2. For Oq ~ [n — {l/2)]/R with large |n|, 
on the other hand, \q'W^ in the twisted field has beats at a frequency about 2A with 
the approximated value of A given in (3.8) in the weak coupling limit, while there is no 
significant beat in the untwisted field with the same parameters, as those in Figure 3. 

3.2 Two-point correlators and detector-field entanglement 

Suppose the initial state of the combined system defined on the t = 0 slice is a factorized 
state, 

IV'(O)) = I 3 ) ® |0l) (3.13) 

which is a product of the ground state of the free detector \g) and the vacuum state of 
the free field |0 l) in the Einstein cylinder x Ri [26]. For the untwisted field, the 
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I 7=0.01, n=2.3.R=21 (-^0.01. £1=2.3, R=21 (-^0.01. £1=2.3, R=2] ( 7=0.01, £1=2,3, R=2 ) 

bif kil’ kSl' I.i1’ 



Figure 5. q'^{t) (upper) and their frequency spectra (lower), n = 3,4, 5,6 (from left to right), 
with the same parameter values in Figure 3 (7 = 0.01, il = 2.3, R = 2). The eigen-frequencies are 
exactly the same as those for the there, but now one can see the beats because there are always 
two dominant eigen-frequencies (w « Oq and |fc„|) for the detector-held mode functions g^"(t). 


vacuum state is further factorized to |0 l) = ® |0l)z) where IOl)^^ is the lowest 

energy state of the field with non-zero wave vector, and \Ql)z is the initial state of the zero 
mode of the free field, chosen to satisfy = 0, which gives the minimal uncertainty 

z( 0 L|(l>fco( 0 ))^| 0 L )2 = z(0L|(nfco(0))^|0L)z = For the twisted field there is no zero¬ 

mode, thus no similar separation is needed. With the factorized initial state (3.13), the 
symmetric two-point correlator of the detectors splits into two parts, e.g., 

{Q\{t)) = {QA{t),QAit)) = I \im{m\QA{t)QA{t') + QA{t')QAit)\m) 

= {QA{t),QA{t))^ +{QA{t),QA{t))^ (3.14) 


where, from (2.13), 


{QA{t),QA{t))^ 

{QA{t),QA{t))y 


h 



E 

ngZ 


h 

2oJri 



2 


(3.15) 

(3.16) 


for the untwisted field, and with kn —)• k'^ for the twisted field. The uncertainty U = 
^J{Q‘a){Pa) ~ {Qa,Pa)‘^ and the purity V = Tr{p^)‘^ = lj{2U), where is the reduced 
state of detector A, will be fully determined by the two-point correlators of the detector 
since the quantum state of the combined system is always Ganssian in this linear system. 


3.2.1 UV cutoff 

In Appendix A, the symmetric two-point correlators of a UD’ detector at rest in R| with 
the initial state similar to (3.13) have been worked out. One can see that {Q\), {Qa^ -^a), 
and {P‘\) all suffer from the UV logarithmic divergence. To control, one has to introduce 
a UV cutoff in the integrals (A.6)-(A.9). The values of cannot be too small, or 
the uncertainty relation of the detector U > h/2 will be violated due to the inconsistency 
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Figure 6 . Time-evolutions of the linear entropy 5'^, contributed by the correlators {P\{t)) and 
{Q\{t)) (upper row) of single detector A in the untwisted field, with the same parameter values as 
those in Figures 2 and 4. Replacing by the gauge-dependent defined in (2.8), both (P^) 
(lower-right) and the corresponding linear entropy S'l^ (lower-left) are always increasing at large time 
scales. The dark-gray, light-gray, and black curves represent the results contributed by the modes 
with wave number = n/R from |n| = 0 up to |n| = 10, 100, and 1000, respectively. One can get 
a taste of the logarithmic divergence from these three curves, with the magnitudes proportional to 
logriniax in the upper-middle plot. The green curves are contributed only by the modes with n = 0 
and \n\ = 3, the latter are almost on resonance with the detector. The blue-dashed and red-dotted 
curves are the green and the black curves with the contributions by the zero mode (n = 0) removed, 
respectively. The difference from other curves with the zero-mode contributions is significant in the 
lower plots. 


with applying the retarded Green’s function in the calculation, similar to the reason we 
mentioned in Section 2. 

How about the correlators for the UD’ detector in x Ri? 

The logarithmic UV divergence from the mode-sum similar to the I 3 term in Appendix 
A will still arise in calculating the two-point correlators of the detector, and will never decay 
out due to the echoes. To get rid of such a UV divergence, again we need to introduce a 
UV cutoff Umax to exclude the modes with |n| > Umax for kn or in our effective theory. 
Again, to be consistent with the retarded Green’s function (2.17) we applied, the UV cutoff 
here cannot be too small to violate the uncertainty relation. 

3.2.2 Linear entropy 

In Figures 6—9 we show selected numerical results of the two-point correlators and the 
detector-field entanglement in this model with the initial state (3.13). The degree of en¬ 
tanglement here is characterized by the linear entropy Sl = 1 — V: the higher the value of 
Sl, the stronger the entanglement between the detector and the field. 

In contrast to the simpler dissipative behavior in Minkowski space, where the time- 
slices are not compact and the spectrum of the field modes is continuous [14, 27], the 
behaviors of the two-point correlators {Q\), {P'^), {QaiPa), and the linear entropy Sl 
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Figure 7. Results similar to those in Figure 6 with the same values of the parameters except 
e = — 1 here for the twisted held. The dark-gray, light-gray, and black curves represent the results 
contributed by the modes with k'^ = [n — (l/ 2 )]/i? from \n\ = 1 up to \n\ = 10 , 100 , and 1000 , 
respectively, while the blue-dashed curve is contributed only by the four modes with |n| = 3 and 4, 
which have the closest frequencies to the detector’s. There seems to be no beat with period longer 
than L in the upper plots. However, the beats with very long period (about 50 x L) emerge in the 
lower-right plot (contributed by the modes with \n\ = 1 to 10). They are produced by the three 
dominant eigen-modes at eigen-frequencies (1.87266,2.30377, 2.75063) [around (wa, Hq, 014 )]. Indeed, 
the beat frequency is about ||2.30377 - 1.87266| - |2.75063 - 2.30377|| « 0.01575 « 27r/(48.84L) 
here. The lower-left and lower-middle plots show the results in terms of the unphysical defined 
in (2.8). The differences from the results with the kinetic momentum are not as significant as those 
in the untwisted field: only the spikes around t = nL in the upper plots are smeared. 


in Si X Ri are more complicated. At early times before the first echo arrives {t/L < 1), 
the correlators and do behave similarly to their Minkowskian counterparts (also see 
Appendix A). At a large time scale, however, the a-parts of the correlators never decay 
out, and the beat of the mode function in Figure 2 is obvious in Figure 6, while there 
is no significant beat in the upper plots of Figure 8, as indicated in Figure 3. This is also 
true for the twisted field in Figures 7 and 9, though the beats of their mode functions are 
not shown in this paper. 

The behavior oi Sl cannot be approximated by including the contributions only by 
the one or two field modes nearest to resonance with the detector in the v-part of the 
correlators. Indeed, the beats in the green and blue-dashed curves in the upper plots of 
Figures 6 to 9 are almost gone, since the contributions of (the upper-left plot in 

Figure 2) and (the second plot from the left in the upper row of Figure 4) have the 

same order of magnitude but are out of phase. Moreover, as we mentioned above, similar 
few-mode approximations can violate the uncertainty relation in some periods of time in 
the system’s history due to the inconsistency with the retarded Green function. During 
these periods, the value of Sl becomes negative and so unphysical. To get rid of this one 
has to include enough field modes in the v-part of the correlators (at least Umax ~ 100 in 
Figures 6—9). 
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The higher-frequency modes are responsible for the small oscillations on top of the beat 
at the frequency about 217 in the weak coupling limit. As the UV cutoff Umax increases, 
such small oscillations of (Q^) and will be amplified, while the whole evolution curves 
of (Q^), and thus Sl, will be elevated as well. Including the higher-frequency modes 

further helps to resolve the spikes of the evolution curve of occurring around t = nL, 
n = 1, 2, 3,... as shown in the middle plots in the upper rows of Figures 6—9. These spikes 
are due to constructive interferences occurring periodically at t = nL when = 1 for 

all = \k^,\ or \k'^i\, n' = 1,2,3,... such that the mode sum of (T*^) get an additional 
logarithmic divergence. More explicitly, the terms in {P\{t))^ corresponding to C 3 (t) in 
(A.5) goes to the counterpart of 13 in x Ri as t —)• nL (in this only occurs as 
t —)• 0). While the value of (F 4 ) varies significantly at these moments, once the UV cutoff 
is introduced, the amplitudes of the spikes relative to their neighborhoods will be finite. 
They can be higher or lower than their neighborhood (see the factor of the C' 3 -term in 
(A. 5)), but will never overwhelm those which have been corresponding to the terms for 
t 7 ^ nL to make the corrections to (^( 4 ) from higher-frequency modes negative. 

In the lower-right plot of Figure 7, while the frequency of the detector is near resonant 
to none of the field mode, the curves happen to show a beat with a very long period. 
Such a beat is produced mainly by three eigen-frequencies rather than two. When the 
free detector frequency is almost, but not exactly, located at the middle point between 
two frequencies of the free field modes, and L has a proper value, the three eigen-modes 
with eigen-frequencies around these three natural frequencies will dominate. While the 
difference between every two dominant eigen-frequencies is 0(1/72), which gives no beat 
beyond the period L of the echoes, the small difference between 1) the frequency difference 
of the middle and the lower modes, and 2 ) the difference of the higher and the middle 
modes, stands out and sets a beat frequency. Note that if the spacing 1/72 in the frequency 
spectrum of the free field is too large, only the eigen-mode with frequency around the free 
detector’s will dominate and so no significant beats will be observed; If 1/72 is too small, 
there will be so many dominant eigen-modes that the evolution becomes complicated. 

3.2.3 Effects of the zero-mode 

For the untwisted field, the curves with or without the zero-mode contribution (e.g. the 
black and red-dotted curves, respectively, in Figures 6 and 8 ) are almost indistinguishable 
in the results with the kinetic momentum Pa- The zero-mode contributions to the v-parts 
of the two-point correlators with Pa are not important here since we have chosen the initial 
state of the zero-mode having the minimal uncertainty. Just like other field modes, a slight 
change of the initial state of the zero-mode will change the results of the correlators and 
Sl slightly. 

When Sj and are replaced by the gauge-dependent S'j and P 4 , the zero-mode will 
make {P'a) growing as indefinitely and so Sl also growing in time at large time scales (see 
the lower plots of Figures 6 and 8 ). This behavior, taken on face value, may be construed as 
the detector-field entanglement increases forever. However, since is a gauge-dependent 
quantity, this ill-behavior is not physically meaningful. Only the entanglement in terms of 
the physical Py^ matters. 
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Figure 8. Another example in the normal (untwisted) field, for comparison with Figures 3 and 
5. The upper plots are the results with the physical and the lower plots are those with the 
unphysical P'j^. The dark-gray, light-gray, and black solid curves represent the results contributed 
by the modes u)k„ from |n| = 0 up to \n\ = 10, 100, and 1000, respectively, while the blue-dashed 
curve is contributed only by the zero mode and the modes with |n| = 4 and 5, which are the modes 
with the frequencies closest to the detector’s. The blue-dashed (|n| = 4,5) and red-dotted curves 
(|n| = 1 to 1000) are obtained from the green and the black curves with the contributions by the 
zero mode removed, respectively. As in Figure 3, no significant beat can be observed here. 



Figure 9. The same parameters as those in Figure 8 except here the detector is in the twisted 
field (e = —1). (Upper row) The dark-gray, light-gray, and black solid curves represent the results 
contributed by the modes = |n— (l/2)|/i? from n = 1 up to n = 10, 100, and 1000 respectively, 
while the blue-dashed curve is contributed only by the modes with n = 5, whose mode frequency is 
the closest to the detector’s. In contrast to the case with the untwisted field in Figure 8, here the 
beating behavior is obvious. (Lower row) (P^) and S'j^ contributed by the unphysical momentum 
P^. The red curves represent (P^) and (the black curves in the upper row) for comparison. 
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Figure 10. Our nonperturbative result for the probability of finding the detector in the first 
excited state (black), is compared with (3.17) obtained by first-order TDPT with finite duration 
of interaction (red dotted) in the perturbative regime (7 <C O). One can see that the deviation 
becomes significant as early as t « L. For comparison, we also show the first-order TDPT result 
(3.18) with infinite duration of interaction weighted by a Gaussian switching function with width t 
(blue). Note that L 7 <C 1 even for t > L in this plot. 


3.3 Comparison with perturbative results 

3.3.1 On the zero mode 

The transition probability of a single UD’ detector in the untwisted field from the initial 
ground state to the first excited state of the detector, obtained using the time-dependent 
perturbation theory (TDPT), depend on the initial two-point correlator of the conjugate 
momentum of the zero mode, (n|^( 0 )) = 2(0i|n|^(0)|0i)2, but is independent of the initial 
correlator of the zero mode itself, (^1^(0)) = z( 0 l|$|^( 0 )| 0 l )2 [24]. Accordingly the authors 
of [24] claimed that the effect of the zero mode in the UD’ detector theory defined in x Ri 
can be suppressed by choosing a suitable initial state such that (11^^ (0)) is very small. They 
pointed out that this is not possible for the usual UD detector with coupling in x Ri 
since in that case ( 0 )) also enters the response function and will become large if (R^q ( 0 )) 
is squeezed (due to the uncertainty relation (4>|^(0))(ft|^(0)) > (fi/ 2 )^). 

From Eqs.(20)-(23) in [27], we have seen that after taking the weak coupling limit 
7 <C D for the complete expression of the transition probability, the leading term will be 
proportional to the 0 ( 7 ) part of {P'^)+Qq{Q\) (mo = 1 here). If we choose the static initial 
state so that z( 0 L|{$fcQ( 0 ), ^[^^(OjjlOL)^ = 0 , the zero mode contributions to {Q\) and (P^) 
will be (Reg^°) 2 ($ 2 ^( 0 )) + and (Rep^°) 2 ($ 2 ^( 0 )) + (Imp|’) 2 (f[|^( 0 )), 

respectively. When 0 < t < L, from (3.5) one sees that both and = d^q^ are purely 
imaginary, which implies that Reg^° = Rep^° = 0 for all time according to Eq.(3.2). Thus 
indeed, the initial variance (^^^( 0 )) of the zero mode are totally irrelevant to the two- 
point correlators of detector A in the untwisted field. The zero-mode contribution to the 
dynamics of the combined system in our model depends on (n|^( 0 )) only. 

3.3.2 On the validity of perturbative transition probability 

In Figure 10, we show the transition probability from the ground state \g) to the first 
excited state |ei) of a UD’ detector initially in the vacuum state of an untwisted field 
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obtained by our nonperturbative method (black), and the first-order TDPT result (red 
dotted) [27] given by 


Pi,i 




(ei|Q^(0)|5) ^ I dh I 


27 I 1 — cos Qot 


0(1^2 


Q2 

“0 


+ 47r ^ 


n 


n=l 


1 — cos -I- Qq) t 


(3.17) 


where h= 1 ^, l>™ (t) 




is the free held operator, and 


1(6110^(0)15)1^ = 1/(200). One can see that the TDPT result (3.17) deviates signihcantly 
(still 0(7)) from the nonperturbative result as early as t ~ L when the hrst echo returns, 
though the echoes appear to be higher-order effects in the equation of motion (3.1). For 
further comparison, we also present the perturbative result (the blue curve in Figure 10) 


Pt 


X 


20I 

7b 


2 roo roo 

dti / dt2Xt(ii)Xi(i2)e"*^‘’(*i"‘2)5^^Q^^^o|l>£Lo(ti)$£Lo(t2)|0) 


0 J —00 
2 


OoL2 


e *^*^0/(8^) + 47r ^ nexp 


n=l 


„ , 27rn ^ 


(3.18) 


with the Gaussian switching function Xti'^) = 2e that the width t = Xt{'^)dT = 

Jq dr corresponds to the duration of interaction in Eq. (3.17). We hnd that the Gaussian 
switching function greatly suppresses the perturbative transition probability (3.18), which 
decays very quickly as the effective interaction time t increases. Gomparing (3.18) with 
(3.17), one can see that at a fixed t, the summand of the summation Yln=T like ne~"'^ 
for large n in (3.18), and n~^ in (3.17). So (3.18) receives much less contributions from the 
short-wavelength modes than (3.17) does if Umax is large enough. 


4 Two-detector case 

Gonsider two identical detectors A and B at rest with natural frequencies uja = ojb = XIq, 
located a.i x = z\ = Q and x = 4 = Rtt = L/2, respectively, as shown in Figure 1 (left) 
and (middle). Substituting the expansions (2.13) and (2.14) into the Heisenberg equations 
of motion for the operators yields 

{df + 2xdt + nl) (t) = 

\ 2 ^ 

n'=l 

(e"' + £“'+') 9 (t - [n' - (1/2)] L) (t - [n' - (1/2)] L)} , (4.1) 

where 5 = A,B,kn (untwisted) or k'^ (twisted field), n G Z, (j}^4 = 41^ = 0 , (44{t) = 
^-i(\k\t-kx) for /c ^ 0 {k can be kn or k'^X), and (44° {t) = 1 — the zero-mode. 

®Note that it should be h~^ instead of h~^ in the overall factor of the RHS of Eq.(26) in [27]. 
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From the setups in the extended coordinates (Figure 1 (middle) and (right)), one 
can see that the classical dynamics for an individual detector and the untwisted field in 
the cases with two or more identical detectors would be equivalent to the dynamics of 
the detector and the untwisted field in the single detector case: in the two-detector case 
QA{t) is affected by Qsit — (-^/2)) with identical solutions to QA{t — {L/2)) in the one- 
detector case. However, this is not true for the mode-functions in quantum theory. In 
(4.1), q^{t) is affected by q^{t) rather than q^{t), and the initial condition for the former 
is q^{t < {L/2)) = 0, which is different from (3.3) for q^{t < {L/2)) = q^{t < (L/2)). 

Suppose the initial state of the combined system at t = 0 is a factorized state, 


V'(O)) = I qA^Qs) ® I Ol), 


(4.2) 


which is a direct product of a two-mode squeezed state of the detectors | <7^, ), or in the 

Wigner function representation, 


Pab{^) — ^ ‘^{Pa + Pb)‘^+ 

0!~‘^{Qa - Qb)'^ + cx^iPA - Pb)'^] 


(4.3) 


with constant parameters a and (3 [14], and the vacuum state of the field | 0^). Again, for 
the untwisted field, we further factorize the field state into | Ol) = | I 0 )z where 

I 0 is the initial state of the zero mode giving z{^L\9k^{^)\^L)z — A^L\^‘kQ{^)\^L)z — 

4.1 Untwisted field 

In the untwisted field, writing q^ = {q^ ± q^)/^/2, the equations of motion (4.1) for the 
mode functions can be simplified to 


{dt + 27^4 -|- Hg) Q'±(^) 


71 ^ 


= -^d, 


± {pit' 9 dtq^^ 

n'=l 




A 






where 


[ 0 ], 


P±{t) = dt (t>Q^{t)±(l)^jt^{t) 




±(/>, 


^ L 


+ 




^ L 


(4.4) 

(4.5) 

(4.6) 


When 0 < t < L/2, F'Tit) = [1 ± (-1)"] for K T 0, Pl°{t) = and 

= F^ = fT = 0. When t > L/2, all the components of F^ vanish. Compared 
with (3.2), one can see that the equations for and q_ in (4.5) are almost equivalent 
to the ones for a single detector in the untwisted and twisted fields, respectively, in 5^ 
with circumference L/2, except that the force terms F^{t) and F^{t) driven by vacuum 
fluctuations of the field in (4.5) are different from their counterpart in (3.2). 

The above delayed differential equations will be solved together with the conditions 
that for 0 < t < L/2, q^{t) = q^it) = 0, g^(t) = q^{t) are given by Eq.(3.3), and 
q\{t) = q^{t) are given by (3.4) and (3.5), before the mutual influences set in. 
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4.1.1 Eigen-frequencies 

Substituting the ansatz ~ , d = A,B into (4.5), we find 

— -|-JIq = — 270 ;+cot(a;+L/4), (4.7) 

—+ 00 = 270 ;-tan(a;_L/4), (4.8) 


which are almost the same as (3.7) and (3.11) in the one-detector case except L there 
has been replaced by L/2 here. The solutions of the above two equations are the eigen- 
frequencies of the mode functions q^. 

In weak coupling limit 7 <C Ho) if is not very close to 2mxjL = n/R, n G N, the 
two detectors with the same natural frequency Oq will be mixed together by the effective 
coupling mediated by the field, and so the dominant eigen-modes will split into 0 ;+ ~ 
Ho -b 7 cot(no-b/ 4 ) from (4.7) and a;_ Hq — 7 tan(noT/ 4 ) from (4.8). These two eigen- 
frequencies can be quite close to each other, and so the beat frequency of = (<?$ ± 

q'^)/V^ will be approximately 


A 


7 



-b tan 




27 

sin(no-b/2)| 


(4.9) 


When rio ~ (2n -b l) 7 r/L, the above beat frequency will reach the lowest value A ~ 27 in 
the weak coupling limit. 

When the natural frequency of the detectors is very close to the frequency of a free 
field mode, namely, Hq ~ 2 n 7 r/L, with some positive integer n large enough, the two 
detectors and the field mode will split into three eigen-modes at frequencies 2 n 7 r/L and 
( 2 n 7 r/L± A), where A ~ will be the beat frequency in this case. This is a half of the 

beat frequency 2 A we obtained below (3.8) in the one-detector case. If the eigen-frequency 
around 2 n 7 r/L is not exactly at the averaged value of the two eigen-frequencies around 
( 2 n 7 r/L ± A), when L has a proper value, there will be beats with an even longer time 
scale similar to the one in Figure 7. In other parameter range of L, the beating behavior 
may disappear or get more complicated. 


4.1.2 Super-radiant mode and instability? 

In Minkowski space, the modes for two detectors similar to here (Eq.(20) in [14]) can 
be interpreted as the super-radiant modes since the dissipation rate is about twice of the 
one for the single detector there. However here in a similar limit, L <C 7 'C Hq, expanding 
q+{t - {L/2)) = q^{t) - {L/2)d^qj^{t) -b ..., Eq. (4.4) gives 


L 

2 



24 J 



^27 -b d/q_^_ -b -b d^q_^_ 


0 , 


(4.10) 


by assuming q_^_ varies slowly in time. This implies d^q+ ~ g(-r±iaj)t p 2 ^ 

HoT^/ 8 )^ <C = Hq -b ( 87 /L). io is nothing but the lowest eigen-frequency obtained from 
(4.7) in the same limit (where cot(Lw/4) ~ 4/(La;)). So q_^_ here oscillates in the lowest 
eigen-frequency rather than decays in a quicker rate. The behaviors of such solutions 
are quite different from the super-radiant modes in [14] since the space here is extremely 
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Figure 11. Entanglement dynamics of two UD’ detectors contributed by the field-modes with kn, 
\n\ = 0,..., 250, in terms of detector ^’s clock. Here 7 = 0.005, H = 2.3, L = Itt, and a and /3 are 
shown in the plot labels. (Upper left) \q^\ (black) has a beating behavior like (l-|-costA)/2 (green- 
dashed) and |q ;4 1^ (gray) like (1 —costA)/2 (blue-dashed) with A « 0.01014 « 27 / sin(UoU/2) given 
in (4.9). (Upper-right and lower row) The logarithmic negativity Ejq- = maxjO, — log 2 2c_} (black) 
of the EnLC compared with the Ejq- of the EnS (gray dashed), the linear entropy S'l of detector 
A (red), and the linear entropies ' of the detector pair evaluated around the future light cones 
(blue) and the hypersurfaces of simultaneity (green dashed) of detector A. The gray part of the 
curve for the £V of the EnLC represents — log 2 2c_ at negative values, where Ej^ = 0. In the lower 
left plot, both of the two detectors are initially in their ground states. One can see that E^ has a 
beat frequency 2A, so does the linear entropy Sl of detector A in the right plots. 


compact: the size of the space L is the same order of the separation of the detector L/2, 
both are extremely small. 

For two UD detectors at rest in (3-1-1) dimensional Minkowski space, when the separa¬ 
tion is small enough, the detectors will become unstable [14]. Now in an Einstein cylinder, 
when the value of L is small, would we get runaway solutions for the mode functions of the 
detector(s)? 

The instability in the analysis in [14] is due to the fact that the amplitude of the 
retarded fields in (3-|-l)D Minkowski space diverges as one approaches to the point-source. 
However, in our (1-|-1)D Einstein cylinder the retarded solution of a massless scalar field 
is regular around the source. Thus a small value of L in our model will produce no 
instability. Indeed, in our calculation we always obtain real eigen-frequencies. Imaginary 
eigen-frequencies never show up when we take L to 0-|-. 
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4.1.3 Entanglement dynamics 

Some examples of the entanglement dynamics for identical detectors A and B located at 
X = 0 and L/2, respectively, in the untwisted field are illustrated in Figure 11. The degree 
of entanglement between the two detectors is characterized by the logarithmic negativity 
[28, 29] = max{0, — log 2 2c_} evaluated around the future light cone (EnLC) of detector 

A, which is essential in quantum teleportation from detector A to detector B [3]: the value 
of the Ej\f of the EnLC is a monotonic function of the upper bound of the optimal fidelity 
of coherent-state teleportation [30]. Here c_ is the lowest symplectic eigenvalue of the 
partially transposed covariance matrix in the reduced state of the detector pair. The 
higher UV cutoff Umax is introduced, or the larger initial variance (n|g) of the zero mode 
is given, the lower EnLC will be. 

In Figure 11 (upper-left), one can see that behaves like (1-|-costA)/2, and \qA^ 
like (1 — costA)/2 with the beat frequency A given in (4.9). Both \q^^ and \q^^ have the 
beat period 2 tiIIS., which is two times longer than the results in the one-detector case in 
Figure 2. In the other plots of Figure 11 with different choices of the parameters a and (5 for 
the initial two-mode squeezed state of the detector pair, the of the EnLC of detector A 
(black curves) in each case shows the same beating behavior of the mode function squared 
P, so do the linear entropy of detector A. One can see entanglement creation, 
sudden-death and revival of the detector pair in the lower plots from the interplay with 
the field. 

The Ej^f evaluated on the hypersurfaces of simultaneity (constant t-slices) (EnS) of 
detector A (gray dashed curves) has a similar beating behavior, but their extrema lag 
behind the EnLC’s in the clock of detector A, so does the linear entropy Sl oi detector A, 
which is a measure of entanglement between a single detector and the rest of the combined 
system (including the field in addition to the other detector). In the right plots, it appears 
that Sl oi detector A increases as the E^ of the ENLC does. However, such a tendency is 
not clear in the lower-left plot, where detectors A and B are initially in a separable product 
state of their ground states. 

( 2 ) 

In the right plots we further introduce the linear entropy of the detector pair, Sj^' = 
1 — Tr(p^g)^ = 1 — (4\/det V), where is the reduced state of the detector pair 

A and B, and V is the covariance matrix of the detector pair. This is a measure of the 
entanglement between the detector pair and the field, combined as a bipartite system. 
Similar to Ej\f, the value of Sj^ can be taken either on the future light cones (blue) or on 
the hypersurfaces of simultaneity (green dashed curves) of detector A. The extrema of the 
latter also lag behind the ones of the former. We find that either way Sj^' does not follow 
the largest beats of E/^. This indicates that the increased entanglement between detectors 
A and B are mainly influenced by the nonlocal correlations in the field, rather than the 
correlations between the detectors and the field. 

4.2 Twisted field 

In the twisted field, and q^ will not influence each other if detectors A and B are 
located exactly at x = 0 and x = L/2 in S^. This is because in detector j4’s point of 
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view, the retarded field from detector B or it’s image at x = (n' — (1/2))L, n' = 1, 2,3,.. 
in the extended coordinates carries a factor, but the field from x = (—n' + (1/2))L 
carries an factor and cancel the former when e = — 1 (note that -L/2 is in the domain 
(—L/2,L/2] for x in the restricted coordinates, but is not). So the situation of our 

setup for the twisted field looks rather simple: = 0 for all time, and the equations 

of motion (4.1) reduce to Eqs. (3.1) and (3.2) for single detectors with the subscript d = A 
generalized to d = A, B while (/>Q'^"(f) for d = in Eq.(3.1) or (3.2) is replaced by ';^^|/ 2 "(^) 
when d = B. Their eigen-frequencies are thus identical to those in Section 3.1.2 with the 
same parameter values. 


4.2.1 Vanishing v-parts of the cross correlators 

The quantum correlation initially in the twisted field will never be converted to the detector- 
detector entanglement since the v-parts of the cross correlators always vanish here. In¬ 
deed, writing = {n — l/2)/R, n = 1,2,3,..., so that for > 0, and 

k'^, = when k'^, < 0, n' = 0,1,2,— Erom Eq.(3.2), one has q^{t) = qj^"{t) 




[ 0 ]-^/ 

= d^(j) Q "(t). Similarly, since 


and p^{t) = p^^"(t) since = —ioo'^e 

(t) = = —f(—l)”'d^0 Q*^"(f) = one can see that 

9b"(^) = and p‘^"(f) = ip^{t) = These imply 

h 


(7^^(^)7^'B(^'))v = 


S 

h 


Re 


{ 


9a (0 




V -\Re\ 

2uL I 


n>0 

h 

n>0 


Re 


eTit) 
{ 9a" (^) 


9b ""(i')} 


sTit') + 

(9b" (^0 + 9b‘^"( i'))} = 0, 


(4.11) 


where (7^, g), (7^^ g') = {Q, q) or {P,p). 


4.2.2 Entanglement dynamics 

Some examples for the dynamics of the EnLC between the two detectors in our setup in 
the twisted field are shown in Eigure 12 with the same values of the parameters as those 
in Eigure 11 except e = — 1. One can see that the entanglement dynamics in this setup 
have the same beating behavior as Iq'aP rather than the combination of |(?aP and I^/aP 
the untwisted field. The beat frequency is the frequency difference between the dominant 
eigen-modes of \q^\'^{= |9bP) around fl, which can be obtained from (3.11). 

The linear entropy Sl of detector A in the upper-right plot of Eigure 12 shows the 
same tendency of increase or decrease as Ej\f of the EnLC, very much like the case with 
the untwisted field (the upper-right of Eigure 11). However, in the lower row of Eigure 12 
the tendency is opposite: Sl increases as the value of Ej\f or — log 2 2c_ decreases. If we 
vary the value of a and /3 continuously, there will be a continuous crossover between the 
former and the latter behaviors, where the Sl curvew have valleys both around the peaks 
and the valleys of E/^. 
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Figure 12. Entanglement dynamics of two UD’ detectors in the twisted field, contributed by 
the field-modes with |n| = 1,..., 1000. All the parameters have the same values as those 
in Figure 11 except e = —1. (Upper-Left) The evolution of |( 7 ^| . {\pa\ ~ \<1a\ this 

weak coupling.) The beat frequency is about 0.0747, which is the frequency difference of the two 
dominant eigen-modes obtained from (3.11). (Upper-right and lower row) The curves represent the 
quantities similar to those in Figure 11. One can see that Ej^/ of the EnLC (black) has the same 
beating behavior as the mode functions. Note that the entanglement creation in the lower-left plot 
for the case with the initial state as the ground states of the detectors are not reliable since the 
linear entropy Sj^ of detector A (red) is negative whenever the Ej^ of the EnLC (black) is positive. 


Unlike the case with the untwisted field in Figure 11, however, here the EnLC and 
( 2 ) 

the S)^ around the future light cones of detector A have opposite behaviors, so do their 

counterparts evaluated on the hypersurfaces of simultaneity. (Again, the extrema of the 

( 2 ) 

Ej\f of the EnS and the 5^ on the hypersurfaces of simultaneity lag behind those evaluated 
around the future light cones in detector A’s clock.) This shows a trade-off between the 
detector-detector entanglement and the field-detector-pair entanglement, due to the lack of 
converting the initial nonlocal field-field correlations to the detector-detector correlations 
in our setup, as we showed earlier in Section 4.2.1. 

Note that the entanglement creation in the lower-left plot in the case with the initial 
state as a separable product of the ground states of the detectors are not reliable here since 
the linear entropy 5^ of detector A is negative whenever the E^r of the EnLC become 
positive. By extrapolation of Umax from our data, we expect that these EnLC will vanish 
if we include enough field modes to make Sj^ of detector A positive at all times (for Umax ~ 
O(IO^), one may have 5^ > 0 almost at all times except some moments around t = nL when 
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the resonance occurs). This would be another difference from the cases with the untwisted 
field as in the lower-left plot of Figure 11, where the detector-detector entanglement can 
be created by the interplay with the field even if the detector pair was started with the 
same initial states. 

5 Summary and Discussion 

5.1 Twisted field, untwisted field, and zero mode 

The Einstein cylinder, with the topology x Ri, possesses two inequivalent configuration 
spaces for a real scalar field [8, 10]. The normal (untwisted) field satisfies the periodic 

boundary condition ^> 3 ;(t) = ^x+L{t) where L is the circumference of S^, and the twisted 
field satisfies the anti-periodic boundary condition <^ 2 ,( 1 ) = —^x+L{t)- 

The untwisted massless scalar field contains a zero mode (with wq = |A:o| = 0), which 
is constant in space [24]. It cannot be simply excluded, otherwise the retarded Green’s 
function of the field will violate causality. Due to its singular normalization, however, 
one has to deal with the zero mode separately from other modes. In terms of the kinetic 
momenta = dtQd associated with the interaction action Sj in (2.3), the effect of the 
zero mode on the dynamics is similar to those from other modes. If we choose Sj in (2.7) 
instead, which yields the canonical momenta + , the zero mode turns out to 

make the two-point correlator of the field amplitudes and the correlators of the canonical 
momenta of the UD’ detectors grow indefinitely, while the energy of the combined detector- 
field system remains constant, and the energy of each subsystem remains bounded. The ill 
behavior in the theory with Sj is simply an illusion when describing the system in terms 
of the gauge-dependent variables P^. 

The dynamics of the detectors in the twisted field is simpler. There is no zero mode, 
and the correlators of the detectors never grow indefinitely. Moreover, in the two-detector 
case, the detector pair in the twisted field will not influence each other if one detector is 
located at x and the other is located exactly at {x -|- L/2 mod L) in S^. In this setup 
the detector-detector mode functions are equivalent to those in the one-detector case, 
and the quantum correlation initially in the field will never be converted to the detector- 
detector entanglement since the v-parts of the cross correlators of the detectors always 
vanish here. 

5.2 Eigen-frequencies and beats 

While the evolution curve of a mode function at early times looks quite similar to the curve 
at late times at a time scale with only few orders of echoes (see Figure 2), as the time scale 
of observation increases, the eigen-frequencies emerge in the frequency spectrum and the 
peaks get sharper. 

The discrete eigen-frequencies of the detector-detector mode functions are the so¬ 
lutions of Eqs. (3.7), (3.11), and Eqs.(4.7-4.8) in the one-detector and two-detector cases, 
respectively. In the weak coupling limit, the eigen-frequencies are close to but not exactly 
the same as the natural frequencies of the free detectors and the free field modes. The 
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eigen-frequencies of the detector-field mode functions are exactly the same as 
though the fluctuations of the field working as driving forces in the equations of motion 
have various frequencies. 

In the one-detector case, when the natural frequency of the free detector flo comes 
close to the frequency of a free field mode, they mix and split into two dominant eigen- 
frequencies, a phenomenon similar to the anti-crossing of energy levels in atomic systems. 
These two neighboring eigen-modes produce the beats in the evolution of a mode function. 
In the weak coupling limit, the beat frequency is determined by the separation of these two 
neighboring eigen-frequencies and proportional to the square root of the coupling strength 
7 . The beats enjoy the largest time scale 0 ( 7 “^/^) in the evolution of the mode functions 
in our compact space S^, which is quite different from the largest time scale 0 ( 7 “^) for 
the detectors in Minkowski space. 

When Qq is not close to any frequency of the field modes, there may not be significant 
beats. However, if Hq is located around, but not exactly at the middle point of two free 
field mode in the spectrum, these frequencies may mix and split into three dominant eigen- 
modes, which may produce the beats at a very large time scale. When this happens the 
beat frequency will be about the difference of the frequency difference between the higher 
and middle eigen-modes, and the frequency difference between the middle and the lower 
eigen-modes. 

In the cases with two identical UD’ detectors, there will be at least two dominant 
eigen-frequencies and so the mode functions always have significant beats, even when Hq is 
not close to any frequency of the field modes. The beat frequency ranges from to 
in the weak coupling limit. 

The beating behavior of the mode functions in turn affects the behavior of the correla¬ 
tors. In the one-detector case, the detector-field entanglement characterized by the value of 
the linear entropy, which is a function of the two-point correlators for the Gaussian states, 
gets the same beating features. This may be interpreted as periodic recoherence. In the 
two-detector case, entanglement dynamics between the two detectors also exhibit the same 
beating feature as in the mode functions. With a suitable choice of the parameter values 
one can observe a sequence of sudden-death and revival of quantum entanglement with a 
(beat) frequency one (twisted field) or two times (untwisted) of the beat frequency of the 
detector-detector mode functions. The beats at large time scale and the resonances around 
t = nL distinguish the dynamics of a detector-field system in x Ri from the dissipative 
behavior of the same system living in Minkowski space R}. 

5.3 Discreteness of field spectrum 

Around the source point, the retarded Green’s function of a field in Minkowski space looks 
the same as the one in any locally flat spacetime of the same dimension. Thus in a classical 
theory, if we couple a source such as a detector with the field in a very short duration, the 
classical physics around the interaction region in two different locally-flat spacetimes will 
be identical. In a quantum theory, however, different field spectra give different Wightman 
functions and produce different vacuum fluctuations. In particular, the field spectrum in a 
non-compact space is continuous, while the one in a compact space or in a cavity is discrete. 
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These differences originating from nontrivial topology or the existence of boundaries can 
affect the dynamics even for a single detector [4, 5, 11]. A comparison of our early-time 
results for a single detector in the Einstein cylinder with the ones in Minkowski space is 
given in Appendix A, where one can see that, while the evolution of the linear entropy 
of a detector in the twisted or untwisted field in the Einstein cylinder looks roughly the 
same as the one in Minkowski space at very early times, in most of the cases studied the 
difference becomes significant as the time scale becomes comparable to the crossing time 
of the compact spatial dimension for the retarded field. In those cases where the natural 
frequency of the detector is in resonance with the frequency of the untwisted held, this 
deviation emerges even earlier, as the evolution curve oscillates about the curve for the 
Minkowski case with a hnite amplitude. 

The discreteness of the held spectrum affects not only the self correlators of a single 
detector, corresponding to the detector-held entanglement, but also the cross correlators, 
which is important in the detector-detector entanglement. Our nonperturbative results 
indicate that early-time evolutions of quantum entanglement between two detectors in 
X Ri and in Rj is small. This is the range where perturbation theory results is expected 
to hold, as shown in Ref. [11] in (3-1-1) dimensions. In our setup the discrepancy between 
the early-time results for the EnLC of a detector pair in 1) the twisted and 2) untwisted 
helds is more pronounced than the discrepancy between the linear entropies Sl of each 
single detector in these two helds. This may be due to the big contrast that the v-parts of 
the cross correlators of the detectors vanish in the twisted held but not in the untwisted 
held. 

We end with a comment on the validity of time-dependent perturbation theory applied 
to entanglement problems in spacetimes with nontrivial topology and boundaries, such 
as studied here. In Minkowski space, TDPT with hnite duration could be valid for the 
interaction time up to 0 ( 1 / 7 ). In x Ri, by counting the exponent of the coupling, we 
can see that the effect of echoes comes from higher-order contributions (e.g., the terms 
in Eq. (3.1)). Nevertheless, our results indicate that even the hrst echo can have signihcant 
effect during L < t < 2L. So we know that as early as t ~ L when the hrst echo returns, 
the lowest order perturbation result may become unreliable (Figure 10). One needs to be 
careful, however, when comparing with the results of [ 11 ] which uses perturbation theory in 
(3-|-l) dimensions with only one spatial dimension being compact. The absence of echoes 
in the results [ 11 ] could be an intrinsic limitation of perturbation theory with a hnite width 
of the Gaussian switching function, or that the echo effect is diluted by the other two non¬ 
compact dimensions. These points are worthy of further investigations by practitioners of 
perturbative switching methods. 
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A Detector in (1+1) D Minkowski space 


For a single UD’ detector at rest in (1-|-1)D Minkowski space, described by the action 
(2.2) and (2.3), initially in its ground state and coupled with the Minkowski vacuum of a 
massless scalar field at t = 0 , the two-point correlators of the detector read 


{Q\t)) = {Q^{t)), + {Q\t)),, 


{P\t)) = {P\t)), + {P\t)), 


with 


{Q\t)), = 


2 UoU 2 

2h'y 
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Figure 13. Early-time linear entropy Sl of a UD’ detector initially in its ground state and 
coupled to the untwisted (blue) and twisted (red dotted) fields in x Ri at f = 0, compared 
with the one in R( (black dashed). Here 7 = 0.01 and rimax = 10000. The natural period of the 
harmonic oscillator is 27r/H « 2.732 in these plots. (Lower right) Frequency dependence of ASj^ = 
I untwisted ~ '5'^ I twisted i which is the difference of Sj^ between the linear entropies in the untwisted 
and the twisted fields at a fixed time t = 0A2L. There are spikes around H = n/R = 2 n 7 r/L, 
n = 1 , 2 ,3,.... 


for ;:g> ri, with n = 0,1, 2,..., T = (0 -|- i'y)t, the sine (cosine) integral function Si 
(Ci), the Euler’s constant 7 e, and the UV cutoff All the above two-point correlators 
of the detector depend on Xg, and so all explicitly depend on the UV cutoff. At late times 
{t ^ 1 / 7 ), one has {Q,P) —)• 0 and {Qp -A 2’h'^X^I'K^ but (P^) —)• 2h'yX^/T: is still cutoff 
dependent significantly. 

The above results are actually identical to the two-point correlators of a UD detector 
in (3-1-1) dimensional Minkowski space [27]. To see this, one may insert = 27rDe^i'’"^® 
to the XgD^ term in while substitute uij^ = 27rDe^“'’"^® to other Xg’s. Thus we 

can borrow the interpretation from Ref. [27] that the constant Ai corresponds to the time- 
resolution of the detector, and the constant Aq corresponds to the time scale of switching-on 
the interaction at the initial moment. 

We have also found that, for the interaction action Sj in (2.7), the correlator {P'‘^) with 
physically non-measurable momentum P' = dtQ + A<ha;=o is both IR and UV divergent. 

In 51 X Pi, the correlators of a single detector at early times before the first echo 
returns can also be expressed as Eqs.(A.2)-(A.5) from (3.3) and (3.4), with the integrals in 
Eqs.(A.6) to (A.9) reduced to the the Riemann sums of the integrands: ^ Yln=T^ 

u ^ ujn = 2n7r/L, n = 1, 2, 3,..., and du —>• Aw = w„+i — Un = 27r/L. In addition, for the 
cases with the untwisted field, one needs to include the zero-mode contribution from (3.5). 
While the difference between each integral of (A.6)-(A.9) in R\ and the corresponding 
discrete sum in x Pi can be large, most of the differences turn out to cancel in the 
correlators. The difference between the correlators with discrete and continuous spectra 


- 29 - 




























can be minimized by fine-tuning the parameters corresponding to the UV cutoff (Aq and 
Ai) in the R\ case. We find that the correlators and Sl in x Ri will not deviate from 
those in R\ significantly until a time scale comparable to L. 

In Figure 13 we demonstrate some examples of the early-time behavior of the linear 
entropy Sl for a UD’ detector initially in its ground state. In each plot, the averaged 
decaying behavior of the three curves at a time scale of 0{L) are very similar, while the 
untwisted field usually gives a greater oscillation at a frequency ss 2VL about a fine-tuned 
curve for the R\ case than the twisted field does. These oscillations are more significant for 
large L (L ^ 27r/fl, upper right), and the most significant on resonance (H ~ 21 x 11 jL with 
integer n, lower row). Fortunately, the amplitudes of such oscillations are always bounded 
and will saturate when t ~ 0 ( 1 / 7 ), provided that t is still not very close to L. To check 
this in the on-resonance cases, in the lower-right plot of Figure 13, we show the frequency 
dependence of the Sl curve with the untwisted field subtracted by the one with the twisted 
field at a fixed time (at each frequency one has some fine-tuned Sl evolution curves in the 
R} case between the untwisted and twisted results). Indeed, one can see that, around the 
resonance peaks, the amplitude of the oscillation at a resonant frequency is not very far 
from those at the neighboring frequencies in value. 

When approaching t = L, in contrast to the curves for the R\ cases, the curves for the 
detectors in the untwisted and twisted fields in x Ri get large but opposite resonant 
behaviors in time, which have the same origin as those spikes around t = nL in Figures 
6 - 12 . 

The evolution curves for the probability of finding the same detector in the first excited 
state, look very similar to the ones in Figure 13. p^i here can be interpreted as the 
transition probability from the initial ground state of the detector to the first excited state 
in TDPT [27]. We find the values of pfi{t) in most of the history in 0 < t < L are about 
0 ( 7 ), indicating that one should be able to see some hints of the non-perturbative p^i{t) 
curves in the leading order (~ 7 ^) of the perturbation theory with a finite duration (cf. 
Figure 10). 
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